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ST ' ABSTRACT. The history of the transition from a neutral intergalactic medium to one that is almost 

^2^' fully ionized can reveal the character of cosmological ionizing sources. In this talk I will discuss the 

' implications for rival reionization scenarios of the rapid decline observed in the space density of quasars 

and star-forming galaxies at redshifts 2^3. The hydrogen component in a highly inhomogeneous 
universe is completely reionized when the number of ionizing photons emitted in one recombination 
time equals the mean number of hydrogen atoms. At z ~ 5, the local character of the UV metagalactic 
flux allows one to define a critical emission rate of hydrogen-ionizing photons per unit comoving volume, 
J^ion = lO^^'^^"'^ phot s~^ Mpc~^. Models based on photoionization by bright QSOs and/or young 
1^ ' galaxies with star formation rates in excess of 0.3 — lM0yr~^ appear to fail to provide the required 

, number of hydrogen-ionizing photons at these redshifts by large factors. If stellar sources are responsible 

OO ■ for keeping the universe ionized at z ^ 5, the rate of star formation per unit comoving volume at this 

i epoch must be comparable or greater than observed at 2: 3. 
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1. Introduction 
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\ The existence of a filamentary, low-density intergalactic medium (IGM) which contains 

■ the bulk of the hydrogen and helium in the universe is predicted as a product of primor- 
dial nucleosynthesis and of hierarchical models of gravitational instability with "cold 

Q ■ dark matter" (CDM) (Cen et al. 1994; Zhang et al. 1995; Hernquist et al. 1996). The 

| H \ application of the Gunn-Peterson constraint on the amount of smoothly distributed neu- 

c/3 ■ tral material along the line of sight to distant objects requires the hydrogen component 

of the diffuse IGM to have been highly ionized by z « 5 (Schneider et al. 1991), and the 
helium component by z « 2.5 (Davidsen et al. 1996). It thus appears that substantial 
sources of ultraviolet photons were present at z > 5, perhaps low- luminosity quasars or 

■ a first generation of stars in virialized dark matter halos with Tvir ^ 10^ K (Couchman 
C3 ; & Rees 1986; Ostriker & Gnedin 1996; Haiman & Loeb 1997; Miralda-Escude & Rees 

1997). Early star formation provides a possible explanation for the widespread existence 
of heavy elements in the IGM (Cowie et al. 1995), while reionization by QSOs may pro- 
duce a detectable signal in the radio extragalactic background at meter wavelengths 
(Madau et al. 1997). Establishing the character of cosmological ionizing sources is an 
efficient way to constrain competing models for structure formation in the universe, 
and to study the collapse and cooling of small mass objects at early epochs. While 
the nature, spectrum, and intensity of the background UV flux which is responsible for 
maintaining the intergalactic gas and the Lya clouds in a highly ionized state at z < 3 
has been the subject of much debate in the last decade, it is only in the past few years 
that new observations have provided reliable information on the presence and physical 



properties of the sources and sinks (due to continuum opacities) of UV radiation in the 
interval 3 < z < 5. 

In this talk I will focus on the candidate sources of photoionization at early times 
and on the timo-dcpondcnt rcionization problem, i.e. on the history of the transition 
from a neutral IGM to one that is almost fully ionized. The starting point of this 
study can be found perhaps in the simple realization that the breakthrough epoch (when 
all radiation sources can see each other in the Lyman continuum) occurs much later 
in the universe than the overlap epoch (when individual ionized zones become simply 
connected and every point in space is exposed to ionizing radiation), and that at high 
redshifts the ionization equilibrium is actually determined by the instantaneous UV 
production rate. In the following I will adopt an Einstein-de Sitter universe (^o = 0.5) 
with Ho = 50/i5o kms~^ Mpc~^ 



2. The Universe After Complete Overlapping 

The complete rcionization of the universe manifests itself in the absence of a Gunn- 
Peterson absorption trough in galaxies and quasars at 2; < 5. In the presence of a 
uniform medium of H i density nui{z) along the path to a distant object, the optical 
depth associated with resonant Lya scattering at Aq;(1 + z) is 

^ -8.3x10 h,, (^^^)3/2' (1) 

where H is the Hubble constant, fa is the oscillator strength of the transition, and all 
other symbols have their usual meaning. The strongest limit on the amount of diffuse 
intergalactic neutral hydrogen at z = 4.3 is provided by a direct estimate of the quasar 
Q1202-0725 flux in regions of the spectrum where lines are absent (Giallongo et al. 
1994): the method yields tqp < 0.02 ± 0.03. Even if 99% of all the cosmic baryons 
fragment at these epochs into structures that can be identified with discrete absorption 
lines, with only 1% remaining in a smoothly distributed component (cf. Zhang et al. 
1998), these measurements imply a diffuse IGM which is ionized to better than 1 part 
in 10^. 

It is useful to start the discussion with the last stage of the rcionization process, when 
individual ionized zones have overlapped, the rcionization of hydrogen (and helium) in 
the universe has been completed, and the IGM has been exposed everywhere to Lyman 
continuum photons. An often overlooked point is worth remarking here, namely the fact 
that even if the bulk of the baryons in the universe are fairly well ionized at 2; < 5, the 
neutral hydrogen still present in the numerous Lya forest clouds and the rarer Lyman- 
limit systems (LLS) along the line of sight significantly attenuates the ionizing flux 
from cosmologically distant sources. In other words, while the complete overlapping 
of I- fronts occurs at Zov > 5, the universe does not become optically thin to Lyman 
continuum photons until much later, at a "breakthrough" redshift Zbr ~ 1-6. It is only 
after breakthrough that any point in space at any redshift z < z^r will be exposed to 
ionizing photons from all radiation sources out to z^t- 



2.1. Cosmological Radiative Transfer 



The radiative transfer equation in cosmology describes the evolution in time of the 
specific intensity J of a diffuse radiation field: 

d a d \ ^ „a ^ ^ c 
ot aou J a 4tt 

where a is the cosmological scale; parameter, c the speed of the light, k the continuum 
absorption coefficient per unit length along the line of sight, and e is the proper space- 
averaged volume emissivity. The mean (averaged over all lines of sight) specific intensity 
of the radiation background at the observed frequency Uo, as seen by an observer at 
redshift Zo, is then 

= i;£</4^i±f (3) 

where v = 1^0(1 + + Zo), and dl/dz is the line element in a Friedmann cosmology. 
The effective optical depth Tes due to discrete absorption systems is defined, for Poisson- 
distributed clouds, as 

re.(.. .) = j^^ dz' I c^iVH: ^^(1 - e-n (4) 

(Paresce et al. 1980) , where d'^N/ dNuidz' is the redshift and column density distribution 
of absorbers along the line of sight, and r is the Lyman continuum optical depth through 
an individual cloud. 



2.2. Intervening Absorption 

The actual distribution of intervening clouds is still quite uncertain, especially in the 
range 10^^'^ < A^hi ^ lO^^'^cm"^, where most of the contribution to the effective 
photoelectric opacity actually occurs. As a function of H i column, a single power-law 
with slope —1.5 appears to provide a surprisingly good description over nearly 10 decades 
in A^Hi (e.g. Hu et al. 1995). At high redshifts, it is a good approximation to use for the 
distribution of absorbers: 

= NoN^,'-'{l + zr. (5) 



dNmdz 



For simplicity, I will assume here a single redshift exponent, 7 = 2, for the entire range 
in column densities. A normalization value of A^'o = 4.0 x 10^ produces about 3 LLS per 
unit redshift at z = 3, as observed by Stengler-Larrea et al. (1995), and, at the same 
epoch, ~ 150 lines above A^hi = 10^^'^''cm~^, in good agreement with the estimates of 
Kim et al. (1997). With this normalization and 7 = 2, the adopted distribution provides 
about the same H i photoelectic opacity as in the model discussed by Haardt & Madau 
(1996). 



2.3. Attenuation Length 

If one extrapolates the N-^-^'^ power-law in equation to very small and large columns, 
the effective optical depth becomes an analytical function of redshift and wavelength, 

4 /-X-1.5 



TeffK,Zo,z) = -V^iVo(^^^ (l + Zo)l-5[(H-z)l-^-(l + Zo)l-^] , (6) 

where uo is the hydrogen photoionization cross-section at the Lyman edge vl, and we 
have not included the contribution of helium to the attenuation along the line of sight. It 
is practical, for the present discussion, to define a redshift ziii'i) such that the effective 
optical depth between Zq and zi is unity. When Vo < J^l, ^ photon emitted at zi with 
frequency i^i — Vo{^ + zi)/{l + Zq) will be redshifted below threshold before being 
significantly attenuated by intervening H i . From equation (^, it can be shown that 
the universe will be optically thin below 

1.616 i/i/i/L 

This expression has a minimum for i^i — 1.59i^l, corresponding to what we shall term 
in the following the "breakthrough epoch" of the universe, Zbr = zi(1.59i^l) = 1.56. 
It is only for z < Zhr that the degree of ionization of the IGM will be determined by 
the balance between radiative recombinations and the total ionizing flux emitted by all 
QSOs which appear after Zbr- 

By contrast, due to the rapid increase in the number of intervening absorbers with 
lookback time, beyond a redshift of 2 the mean free path of photons at 912 A becomes 
so small that the radiation is largely "local" , as sources at higher redshifts are severely 
absorbed. Expanding equation around 2 = 3, for example, one gets Tcg{i'L) ~ 
0.36(1 -I- z)^Az — 1 for Az — 0.18. This corresponds to a proper distance or "absorption 
length" of only 

/ 1 I \ -4.5 

A/K)«33Mpc^ . (8) 



In the local (or "source-function" ) solution to the equation of radiative transfer, this 
strong attenuation effect is approximated by simply ignoring sources with z > zi and 
neglecting absorption for those sources with z < zi.^ Since only emitters within the 
volume defined by an absorption length contribute significantly to the background in- 
tensity shortward of the Lyman edge, cosmological effects such as source evolution and 
frequency shifts can be neglected, and the number of ionizing photons per unit proper 
volume at z can be written as 

An p , Jiiy) . (AO 
nion = — / di^—, — = n-ion , (9) 



^ As filtering through a clumpy IGM will steepen the UV background spectrum (Miralda-Escude 
& Ostriker 1990), the absorption length at the He II edge will be significantly smaller than at 1 ryd. 



where nion(i) = d,i'e{i',t)/hi', and (AZ) is an average over the the incident photon 
spectrum.^ The smaU absorption length is mostly due to systems with continuum optical 
depth r ~ 1. Within the assumption that all absorbers are optically thin, highly ionized 
in both H and He, and contain most of the baryons of the universe, it is possible to 
derive a relation between the mean absorption length and the gas recombination time 
tree from the equation of ionization equilibrium, 

(AO^^cW, (10) 
nil 

where nn is the volume-averaged hydrogen density of the expanding IGM, nii{0) = 

1.3 X 10-^ (nbhl^/om), 

W = [(l + 2x)nHaBq-^ = 1.5Gyr(^^) ^ (^) ^ ' ^^^^ 

as is the recombination coefficient to the excited states of hydrogen, x the helium to 
hydrogen cosmic abundance ratio, C = ('T'hii)/^hii ^1^^ ionized hydrogen clumping 
factor,^ and we assumed a gas temperature of 10^ K. An empirical determination of 
the dumpiness of the IGM at high redshifts is hampered by our poor knowledge of 
the ionizing background intensity and the typical size and geometry of the absorbers. 
Numerical N-body /hydrodynamics simulations of structure formation in the IGM within 
the framework of CDM dominated cosmologies have recently provided a definite picture 
for the origin of intervening absorption systems, one of an interconnected network of 
sheets and filaments, with virialized systems located at their points of intersection. In the 
simulations of Gnedin & Ostriker (1997), for example, the clumping factor rises above 
unity when the collapsed fraction of baryons becomes significant, i.e. z < 20, and grows 
to C > 10 (40) at z « 8 (5) (because of finite resolution effects, numerical simulations 
will actually underestimate clumping). The recombination time is then much shorter 
than that for a uniform IGM. 

3. Reionization of the Universe 

In the rest of this talk I will show how a simple formalism can shed some light on the 
time-dependent reionization process, and assess the role of quasars and star-forming 
galaxies as candidate sources of photoionization at high redshifts. 

3.1. Time Evolution of H ii Filhng Factor 

When an isolated point source of ionizing radiation turns on, the ionized volume initially 
grows in size at a rate fixed by the emission of UV photons, and an ionization front 
separating the H ii and H i regions propagates into the neutral gas. Most photons travel 

2 Since Al{u) oc (u/ul)'^-^ , one has {Al) = Al(i^L)S{'i^-^~"= where 5 = Qs/(1.5-Os), 

e(i^) oc , and we have assumed a cutoff at 4 ryd because of He II absorption. A spectrum with 

as = 1.8 yields (Al) = 2.2AI(ul). 

^ This may be somewhat lower than the total gas clumping factor if higher density regions are less 
ionized (e.g. Gnedin & Ostriker 1997). 



freely in the ionized bubble, and are absorbed in a transition layer. Across the I-front 
the degree of ionization changes sharply on a distance of the order of the mean free path 
of an ionizing photon. When tree t, the growth of the H ii region is slowed down by 
recombinations in the highly inhomogeneous IGM, and its evolution can be decoupled 
from the expansion of the universe. Just like in the static case, the ionized bubble will 
fill its time-varying Stromgren sphere after a few recombination timescales. 

The filling factor of H ii regions in the universe, Qhii, is, at any given instant t, equal 
to the integral over cosmic time of the rate of ionizing photons emitted per hydrogen 
atom and unit cosmological volume by all radiation sources present at earlier epochs, 

/* !^i2Il^dt', (12) 
Jo ^H(t') ^ ^ 

minus the rate of radiative recombinations, 

Qmi{t') 



Differentiating, one gets 



^icc(^ ) 

^Qhii nion Qmi 



dt'. (13) 

(14) 



dt nn iroc 

It is this simple differential equation - and its equivalent for expanding helium zones - 
that statistically describes the transition from a neutral universe to a fully ionized one, 
independently, for a given emissivity, of the complex and possibly short-lived emission 
histories of individual radiation sources, e.g., on whether their comoving space density 
is constant or actually varies with cosmic time (cf Arons & Wingert 1972). In the case 
of a time-independent clumping factor, equation (^4|) has formal solution 

gHii(i) = /* dt' ^ exp (-^ + -^Y (15) 

with <rcc oc t'^ . At high rcdshifts, and for an IGM with C ^ 1, one can expand around 
t to find 

Qmi{t) « —tree (16) 

The porosity of ionized bubbles is approximately given then by the number of ionizing 
photons emitted per hydrogen atom in one recombination time. In other words, because 
of hydrogen recombinations, only a fraction t^cc/t {'^ a few per cent at 2: = 5) of the 
photons emitted above 1 ryd is actually used to ionize new IGM material. The universe 
is completely reionized when Q = 1, i.e. when 

riionticc = n-H- (17) 

Figure 1 shows the H 11 (and He iii ) filling factor as a function of redshift for a QSO 
photoionization model and a clumpy IGM with C = 1,10,20 and 35. The transition 
from a neutral to a ionized universe takes place too late in this model, as late as 2: = 3.8 
for C = 35, and never before z — 4.9 even in the limiting case of a uniform medium. 



redshift 



Fig. 1. The evolution of the H ii {solid lines) and He iii {dashed lines) QWing factors as a function 
of redshift, for a universe where photoionization is dominated by QSOs (from Haardt et al. 
1998). The total IGM density is taken to be ^ibh^o = 0.06, and, from right to left, the four 
curves assume a clumping factor C — 1, 10,20, and 35. The QSO intrinsic spectrum varies as 
* shortward of the hydrogen Lyman edge. Note how the ionization of H i (He ii ) is never 
completed before z = 4.4 {z = 3.5) in models with C > 10. 

3.2. Critical Emissivity 

We can better quantify now the inadequacy of known sources of radiation to provide the 
required amount Lyman continuum photons at z « 5. Because of the local character of 
the UV background at high rcdshifts, one can compute at any given epoch a critical value 
for the photon emission rate per unit cosmological comoving volume, A/ion, independently 
of the (unknown) previous emission history of the universe: only rates above this value 
will provide enough UV photons to keep the IGM ionized at epoch. One can then 
compare our determinations of Mon to the estimated contribution from QSOs and star- 
forming galaxies. We have seen in the previous sections that only photons emitted within 
one recombination timescales or, equivalently, within one attenuation length of any fluid 
element in the universe can actually be used in the equation of ionization equilibrium. 
Equation ( p^ can then be rewritten as rtion(AZ) — tihc, or 



"-h(0)c 



(10^* phot s"^ Mpc"3) /i5o(l + z) 



4.5 ^bh^Q 

0.06 ' 



(18) 



ion — 



2.2AI{i^l) 



where a spectral index of as = 1.8 has been assumed for the volume emissivity between 
1 and 4 ryd. The uncertainty on this critical rate is difficult to estimate, as it depends 
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Fig. 2. Comoving emission rate of hydrogen Lyman continuum photons (solid line) from QSOs, 
compared with the minimum rate (dashed line) which is needed to fully ionize an Einstein-de 
Sitter universe with /150 = 1 and fit — 0.06 (from Haardt et al. 1998). Models based on 
photoionization by quasar sources fall short by about an order of magnitude at z = 5. The 
starred symbols show the estimated contribution of star-forming galaxies at z = 3 and, with 
larger uncertainties, at 2 = 4. They should be considered as upper limits, as they have been 
derived assuming a negligible amount of intrinsic H i absorption. 

on the statistics of intervening absorbers and the nucleosynthesis constrained baryon 
density. A quick exploration of the available parameter space indicates that the uncer- 
tainty on A/ion is unlikely to exceed ±0.3 in the log. The evolution of the critical rate as 
a function of redshift is plotted in Figure 2: A/ion is higher than the quasar contribution 
already at 2: > 3. At z = 5, the deficit of Lyman continuum photons is about an order 
of magnitude. For relatively bright galaxies to produce enough UV radiation at z = 5, 
their space density would have to be comparable to the one observed at z « 3, with 
most ionizing photons being able to escape freely from the regions of star formation 
into the IGM. This scenario appears quite improbable, in light of the decrease in the 
UV galaxy emissivity observed in the HDF above z = 2 (Madau et al. 1996) and of 
direct observations of local starbursts below the Lyman limit showing that at most a 
few percent of the stellar ionizing radiation produced by these sources actually escapes 
into the IGM (Leitherer et al. 1995). 

It is interesting to convert the derived value of A/ion at z = 5 into a star formation 
rate per unit (comoving) volume, p: 

p(z = 5) = Mo„ X 10-^3.i_^^--i ^ 0.025/,-! Mq yr-i Mpc'^, (19) 



where f^sc is the unknown fraction of Lyman continuum photons from O stars which 
escapes the galaxy H i layers into the IGM. The conversion factor assumes a Salpeter 
IMF with solar metallicity, and has been computed using Bruzual & Chariot (1998) 
population synthesis code. If dwarf galaxies associated with cr w 20 km s"-'^ dark matter 
halos are responsible for ionizing the universe at z = 5, and are able to convert a fraction 
fSlb of their total mass into stars over a timescale of order the free-fall time of the halo, 
it can be shown then that their comoving space density must be approximately 

0.025 /-jMpyr-^Mpc-^ ^ 0.13 , 

0.2/Moyr-i W ' ^ ' 

at least several hundred times larger than the space density of present-day L > L* 
galaxies. 
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